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rawinsonde  height  measurements  of  wind  and  temperature  are  fit  using  a 
Hermite  Interpolation  algorithm.  Richardson  numbers  (Ri)  calculated  from 
these  profiles  are  analyzed  to  determine  the  number  and  percent  of 
occurrences  of  Rl  -  Ric  on  a  seasonal  basis.  Summary  (seasonal)  statistics 
on  this  critical  Richardson  number,  or  turbulence  indicator,  as  well  as  all 
variables  used,  are  computed  for  1  km  altitude  bins  from  2  to  25  km,  a 
range  determined  to  provide  good  statistical  accuracy  and  precision.  Par¬ 
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analyses  demonstrate  that  the  patterns  of  turbulence  based  on  a  critical 
Richardson  number  of  1  have  a  substantial  component  which  is  stable  and 
reproducible  from  year  to  year.  The  regressions  relate  percent 
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range  from  0.62  to  0.7  8  (that  is,  up  to  60%  variation  in  percent  occurrences 
explained).  Yearly  variations  increases  the  square  of  the  multiple  correla¬ 
tion  coefficients  only  by  a  maximum  of  0.  003  (that  is,  at  most  3%  variation 
explained). 
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Modelling  Turbulence  in  the  Lower  Atmosphere 
Using  Richardson's  Criterion 


1.  INTRODUCTION 

Elaborate  and  expensive  experimental  techniques  are  frequently  used  to  obtain 
measurements  for  atmospheric  turbulence  studies.  A  common  approach  to  tur¬ 
bulence  investigation  is  to  obtain  accurate  measurements  of  wind  velocity  fluctuations 
for  spectrum  analysis.  There  are,  however,  many  subtle  problems  involved  in  the 
practical  analysis  of  finite  lengths  of  data.  In  the  case  of  finite  lengths  of  periodic 
data  where  the  fundamental  wavelength  is  known  and  digitization  is  performed 

properly,  there  are  no  problems  associated  with  interpretation  of  the  resultant 

12  3 

spectra.  ’  This  is  not  the  case  with  turbulence;  as  pointed  out  by  Hinze,  in  his 

Introduction;  "Turbulence,  consists  of  many  superimposed  quasi-periodic  motions. 

Though  a  harmonic  analysis  of  the  velocity  fluctuations  can  be  carried  out,  this 

fact  is  no  proof  that,  conversely,  the  turbulent  fluctuations  are  composed  of  these 

harmonics.  " 

The  approach  we  have  used  here  consists  of  a  statistical  analysis  performed  on  a  large 
data  base  of  twice  dai  ly  rawinsonde  measurements  of  w  ind  and  temperature  in  the  altitude 

(Received  for  publication  14  December  1981) 

1.  Oppenheim,  A.V.,  and  Schafer.  R.  W.  (1975)  Digital  Signal  Processing, 

Prentice -Hall,  lrc.  ,  Englewood  Cliffs,  New  Jersey.  " 

2.  Brigham,  E.A.  (1974)  The  East  Fourier  Transform,  Prentice-Hall,  Inc., 

Englewood  Cliffs,  New  Jersey. 

3.  Hinze,  J.O.  (1959)  Turbulence,  McGraw-Hill,  Inc.,  New  York,  N.Y. 
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range  2  to  25  km.  The  Richardson  criterion  (Section  2)  is  applied  to  dally  data  and 
the  likelihood  or  relative  frequency  of  occurrence  of  turbulence  baaed  on  a  critical 

Richardson  number  of  one  (Ri  =  1)  ia  determined  on  a  seasonal  basis.  A  similar 

c 

study  of  Richardson  numbers  in  the  mesosphere  using  rocket  grenade  data  has  been 

4 

reported  by  Zimmerman  and  Murphy.  These  data,  however,  were  only  available 
at  six  locations  and  the  number  of  measurements  obtained  were  only  suitable, 
statistically,  for  determining  semiannual  variations.  Rawinsonde  data  used  in  this 
study  have  the  advantage  that  they  are  available  on  magnetic  tape  from  many  stations 
at  low  cost  (140  locations  of  which  90%  are  in  the  Americas).  Also,  data  collection 
techniques  are  standardized  and  all  data  are  edited  at  a  central  location;  the  National 
Climatic  Center.  Seasonal  height  profiles  of  the  occurrences  of  turbulence  based 
on  a  critical  Richardson  number  of  1  have  been  qualitatively  evaluated  at  many  of 
these  stations.  These  height  profiles  revealed  interesting  structure  and  variability 
with  altitude,  latitude,  and  season.  These  features  will  be  presented  in  this  report 
along  with  a  regression  analysis  performed  on  data  from  21  of  the  stations  located 
within  the  continental  United  States.  The  objective  was  to  determine  the  significance 
of  a  relationship  between  the  proportion  of  Richardson  numbers  less  than  unity 
(turbulence  measure),  and  altitude,  location,  and  season. 
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Figure  1.  Occurrences  of  Turbulence,  The  percent  occurrences  are 
obtained  by  taking  the  ratio  of  the  total  number  of  occurrences  of 
Ri  £  Ric  to  the  total  number  of  measurements  at  a  given  altitude.  This 
is  repeated  for  Ri  ranging  in  steps  from  0.  25  to  1.  00 


4.  Zimmerman,  S.  P.  ,  and  Murphy,  E.A.  (  1977)  Stratospheric  and  mesospheric 
turbulence.  Dynamical  and  Chemical  Coupling,  B,  Grandal  and  J.A.  Holtet, 
Eds..  D.  Keitel  Publishing  Co.  ,  Dordrecht,  Holland,  AFGL-TR-78-0020, 
AD  A050128. 
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Figure  2.  Location  of  the  140  Available  Kawinsonde  Reporting  Stations 


Table  1.  Rawinsi  nde  Data  Collection  Stations 


Name 

Latitude 

1 

Longitude 

ATLANTIC  STATION  H 

38 

00N 

071 

OOW 

SAN  DIEGO.  CA 

32 

4  9N 

117 

08VV 

YUCCA  FLAT,  NV 

36 

57N 

116 

03  W 

SAN  NICOLAS.  CA 

33 

16N 

119 

33W 

HUNTINGTON,  WV 

38 

o  2\T 

082 

32\V 

SALEM,  IL 

38 

39N 

088 

58W 

CENT REVIL.LE.  AL 

32 

54N 

87 

1 5  W 

LAKE  CHARLES,  LA 

30 

07  N 

093 

13W 

JACKSON.  MS 

32 

19N 

909 

05  W 

MONETT,  MO 

36 

53N 

093 

54  W 

LONGVIEW.  TX 

32 

2  IN 

094 

39W 

BOGOTA,  COLUMBIA  SA 

04 

4  2N 

074 

90W 

BARBADOS,  \V1 

13 

04  N 

059 

30W 

SANTO  DOMINGO,  DOM.  HEP. 

18 

28N 

069 

53  W 

. JU!iililluiiiiil.i.A 


Table  1.  Rawlnaonde  Data  Collection  Stations  (Cont) 


Name 

Latitude 

Longitude 

TRINIDAD,  WI 

10 

35N 

061 

21 W 

SAN  JUAN,  PR 

IB 

26N 

066 

OOW 

WILLEMSTAD,  CURACAO  N.  A. 

12 

1  IN 

058 

58W 

ST.  MAARTEN.  N.  A. 

18 

03N 

063 

07W 

GUANTANAMO,  CUBA  NAS 

19 

54N 

075 

09W 

KINGSTON,  JAMAICA  W  I 

1? 

56N 

076 

47W 

SWAN  ISLAND 

1? 

24N 

083 

56W 

GRAND  CAYMAN.  WI 

10 

18N 

081 

22W 

SAN  ANDRES,  COLOMBIA 

12 

35N 

081 

4  2W 

CHOLUTECA,  HONDURAS 

13 

18N 

087 

11W 

GUATEMALA  CITY,  GUATEMALA 

14 

32N 

090 

34W 

MEXICO  D  F  1  A 

19 

26N 

099 

04W 

VERACRUZ.  MEXICO 

19 

ION 

096 

C7W 

APALACHICOLA,  FL 

29 

44N 

085 

02W 

MIAMI,  FL 

25 

48N 

080 

16W 

TAMPA  BAY,  FL 

27 

4  2N 

082 

24W 

KEY  WEST,  FL 

24 

35N 

081 

4  1W 

MERIDA  YUCATAN.  MEXICO 

20 

57N 

089 

39W 

BOOTH VILLE.  LA 

29 

20N 

089 

24  W 

victoria,  tx 

28 

5  IN 

096 

55W 

BROWNSVILLE,  TX 

25 

54N 

097 

26  W 

ST.  GEORGE.  BERMUDA 

32 

22N 

064 

41W 

GREENSBORO.  NC 

36 

05N 

07  9 

57  W 

DAYTON.  OH 

39 

52N 

084 

07W 

WAYCROSS,  GA 

31 

15W 

082 

24W 

ATHENS.  GA 

33 

57N 

083 

19W 

CHARLESTON.  SC 

32 

54N 

080 

02W 

NASHVILLE.  TN 

36 

15N 

036 

34W 

STEPHENVILLE,  TX 

32 

13N 

098 

11W 

LITTLE  ROCK,  AR 

34 

50N 

092 

16W 

OKLAHOMA  CITY.  OK 

35 

24N 

097 

36W 

DODGE  CITY.  KS 

37 

46N 

099 

58W 

TOPEKA.  KS 

39 

04N 

095 

37  W 

CARIBOU.  ME 

46 

52N 

068 

01W 

CHATHAM,  MA 

41 

40N 

069 

58W 

BUFFALO.  NY 

42 

56N 

07  8 

44W 

ALBANY.  NY 

42 

45N 

07  3 

48W 

PORTLAND.  ME 

43 

39N 

07  0 

19W 

FLINT.  MI 

42 

58N 

083 

44W 

PEORIA.  IL 

40 

4  ON 

089 

4 1 W 

SAULT  STE  MARIE.  Ml 

46 

28N 

084 

22W 

GREEN  BAY.  WI 

44 

29N 

088 

08W 

INTERNATIONAL  FALLS.  MN 

48 

34N 

093 

23W 

ST.  CLOUD.  MN 

45 

33N 

094 

05  W 

HURON,  SD 

33 

12N 

098 

13W 

KEFLAVIK.  ICELAND  FWF 

63 

58N 

022 

36W 

COCORRO  IS.  ,  MEXICO 

18 

14N 

111 

03W 

HILO.  HI 

19 

43N 

155 

04  W 

JOHNSTON  ISLAND 

16 

44N 

169 

3 1 W 

CHIHUAHUA,  MEXICO 

28 

4  2N 

106 

04  W 

MAZATLAN,  MEXICO 

23 

1  IN 

106 

25W 

DEL  RIO,  TX 

29 

22N 

ICO 

55W 

MONTERREY.  MEXICO 

25 

52N 

100 

12W  j 

EMPALME,  MEXICO 

27 

57N 

110 

48W 

Table  1.  Rawinsonde  Data  Collection  Stations  (Cont) 


Name  Latitude  Longitude 


GUADALUPE  IS.  .  MEXICO 

28 

53N 

1  18 

18W 

L1HUE  KAUAI.  HI 

21 

5PN 

159 

21W 

BAKKLNG  SANDS.  HI 

22 

02N 

159 

47  W 

MIDWAY  ISLAND.  NS 

28 

13N 

177 

2 1 W 

MIDLAND.  TX 

31 

56  N 

102 

1 2W 

EL  PASO.  TX 

31 

49N 

106 

24  W 

AMARILLO.  TX 

35 

14  N 

101 

4  2W 

ALBUQUERQUE,  NM 

35 

03N 

106 

37  W 

DENVER.  CO 

39 

46  N 

104 

5  3  W 

GRAND  JUNCTION,  CO 

39 

07  N 

108 

3  2W 

ELY.  NY 

39 

17  N 

1 14 

5  1 W 

TUCSON,  AZ 

32 

07  N 

110 

56  W 

WINSLOW,  AZ 

35 

0 1 N 

110 

44  W 

OAKLAND,  CA 

3? 

45N 

122 

13W 

UISMARCK.  ND 

46 

46  N 

100 

45  W 

LANDER,  \VY 

42 

49N 

108 

44  W 

NORTKPLATT'E.  ne 

4  1 

03N 

100 

4  1VV 

RAPID  CITY.  SD 

44 

03  N 

103 

04  W 

SALT  LAKE  CITY.  UT 

40 

46N 

1  1 1 

58W 

\\  INN  EMU  CCA.  NV 

40 

54  N 

117 

48W 

UOISE.  ID 

4  3 

34  N 

116 

1 3  W 

GREAT  FALLS.  MT 

4? 

29N 

1 1 1 

2 1 W 

SPOKANE,  W'A 

47 

38N 

117 

32W 

MEDFORD,  OR 

42 

22N 

122 

5  2W 

SALEM.  OR 

44 

55.N 

123 

01W 

ANNETTE,  AK 

55 

02N 

131 

34  W 

YAKUTAT,  AK 

59 

3  IN 

139 

4  0W 

KODLAK,  AK 

57 

45N 

152 

29W 

KING  SALMON.  AK 

58 

4  IN 

156 

39W 

COLD  HAY.  AK 

55 

12N 

162 

43W 

ADAK.  AK 

51 

5  3N 

176 

39W 

ST.  PAULIS..  AK 

57 

09N 

170 

1 3W 

ANCHORAGE,  AK 

61 

10N 

150 

01 W 

FAIRBANKS.  AK 

64 

49N 

147 

52W 

MCGRATH.  AK 

62 

58N 

155 

37W 

BETHEL.  AK 

50 

47N 

161 

4  8W 

KOTZEBUE,  AK 

66 

5  2.N 

162 

38W 

NOME,  AK 

64 

3  ON 

165 

26W 

BARTER  ISLAND.  AK 

70 

08N 

143 

38W 

BARROW.  AK 

7  1 

18N 

156 

47  W 

YAP  CAROLINE  IS. 

09 

29N 

138 

05E 

KOROR  CAROLINE  IS. 

07 

20N 

134 

29E 

PONAPE  CAROLINE  IS. 

06 

58N 

158 

13E 

truk.  Caroline  is. 

07 

28N 

151 

5  IE 

KWAJALEIN  ATOLL  MI  PMP  WS 

08 

44N 

167 

44E 

MAJURO,  MARSHALL  IS. 

07 

05N 

171 

23E 

GUAM  TAGUAC  MI 

13 

33N 

144 

50E 

WAKE  ISLAND 

19 

17N 

166 

39E 

LIMA,  PEHU 

12 

01S 

077 

08  W 

ANTOFAGASTA,  CHILE 

23 

25S 

07  0 

28W 

QUINTERO.  CHILE 

32 

47S 

071 

3  2W 

PUERTO  MONIT.  CHILE  EL  TEMPUAL 

41 

26S 

07  3 

07  W 

PUNT  A  ARENAS,  CHILE 

53 

02S 

070 

51W 

PAGO  PAGO.  AMERICAN  SAMOA 

14 

20S 

170 

43W 

T.-ble  ..  Rav.  insonde  Data  Collection  Stations  (Conti 


Nall*  C 

L  itit 

-ide 

Longi 

tud<_ 

P 

L  \S1LK  bLAM,  J1I.1-: 

27 

10S 

109 

26  U 

m  «U5  ST  A.  A  \  .  RC  .  1C 

f.O 

01S 

;  ;  j 

3 1 U 

DIEGO  GARCIA 

07 

13S 

07  2 

24  E 

MCMURDO  SOUND  ANTARCTIC 

77 

53S 

166 

44E 

1 1ALLET T  A  NT  A  K  U  Ic. 

AAL'NDSEN  SCOTT  ANTARCTIC 

72 

19S 

170 

13E 

CHINA  LAKE,  CA 

35 

4  IN 

117 

41W 

PT  MUGU,  CA 

34 

06  N 

119 

07  U 

SAN  NICOLAS.  CA  PMR  WS  SITE  1 

33 

14N 

119 

27  W 

VANDENBERG  AFB,  CA 

34 

45N 

120 

34  W 

CAPE  HATTERAS,  NC 

35 

16N 

07  5 

33  W 

STERLING,  VA 

38 

59N 

07  7 

28W 

WALLOPS  IS.  ,  VA 

37 

5  1 N 

07  5 

?9W 

GLASGOW.  MT 

48 

13N 

106 

37  W 

QUILLA  V  UTE,  WA 

47 

57.N 

124 

33W 

NEW  YORK.  NY 

40 

47N 

07  3 

4GW 

PITTSBURGH,  PA 

40 

32N 

080 

14W 

OMAHA,  NE 

41 

22N 

096 

01W 

2.  THE  RICHARDSON  CRITERION 


5 

Here  wc  will  make  use  of  the  criteria  for  stability  developed  by  Richardson. 
Consider  his  energy  budget  equation  for  a  large  volume  of  the  atmosphere  (where 
he  has  assumed  that  the  eddy  viscosity  p  =  eddy  conductivity  c ) : 


2  2 


r  /  0  v  \ 

/  0~  \ 

gao] 

1  „  *  + 

(rf) 

dZ  dA  . 


(1) 


The  second  term  or  the  left  side  is  the  dissipation  due  to  viscosity  and  the  terms 
on  the  right,  the  supply  from  the  mean  wind  and  the  loss  from  convective  activity. 
With  the  assumption  that  the  eddy  viscosity  p  equals  the  eddy  conductivity  c, 
Richardson  considered  a  special  application  of  Eq.  (1)  where  the  volume  of  atmos¬ 
phere  is  initially  at  rest  then  set  into  motion  by  an  approaching  depression  so  that 
c  and  E'  are  just  somewhat  greater  than  zero.  Under  these  conditions  with  the 
dissipation  term  negligibly  small,  then  E'  will  increase  if 

(i£)  *(■£■) »  f  a.  <*> 

p 


5.  Richardson,  L.  B.  (1920)  The  supply  of  energy  from  and  to  atmospheric  eddies, 
Proc,  Roy,  Soc,  A.  j)7:354-373. 
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where  ^  is  the  entropy  gradient.  - 

From  the  first  law  of  thermodynamics 

dq  =  dE  ♦  PdV  .  _ _ (3) 


Since 


and  the  internal  energy  for  an  ideal  gas  is  7ero,  then 

dE  ■(!?-)  dT  =  cvdT 
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(5) 
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and 

dq  =  Cy  dT  *  PdV  .  (6) 

Differentiating  the  f  as  equation 

PdV  <  VdP  *  RdT  (7) 

and  substituting  into  Eq.  (6)  for  PdV,  results  in  the  following  form  of  the  first  law: 

dq  =  (C  +  R)  dT  -  VdP  . 
v 

Since 

C  -  C  =  R 
P  v 

then 

dq  -  C  dT  -  VdP  . 

P 

Substituting 
V  =  RT/P 

and  dividing  by  the  integrating  factor  T  to  obtain  an  exact  differential,  the  change  in 
entropy  is  defined  as 


(8) 


(9) 


(10) 


(11) 
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The  constant  ■&- 


■( 


-  dT 
ST 


can  be  obtained  from  a  differentiated  form 

1  f 

P  '  adiabatic 

of  Poisson's  equation  with  the  potential  temperature,  9  ,  independent  of  altitude. 

,  the  dry  adiabatic  lapse  rate, 


With  the  definition  =  •£-  -  ^  ^  j 


P  '  'adiabatic 
then  tne  original  criteria,  Eq.  (2),  can  be  written  as  a  ratio 
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dT 


♦  r 


ISZ  a  \ 


(5T  *  (£ 
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(15) 


Also,  logarithmic  cifferentiation  of  Poisson's  equation  with  the  potential 
temperature  varying  w"'  ight  yields  the  relationship: 


1  a  9  1  f  dT  gl 

V  FZ  ‘  T  L  ST  +  cj 


(16) 


Except  for  a  factor  g,  this  is  equivalent  to  the  numerator  in  Eq.  (15)  and  yields  the 
other  familiar  form  of  the  criteria: 
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A  form  of  the  stability  criterion  was  first  used  and  named  after  Richardson 
6  T 

by  Paeschke.  Swinbank  was  somewhat  critical  of  Richardson's  earlier  work  on 
the  conditions  for  the  onset  of  turbulence.  In  a  heated  reply  by  Richardson  to 
Swinbank  which  follows  Swinbank's  paper,  Richardson  mentions  his  embarrassment 
at  having  to  note  Paeschke' s  failure  to  take  Into  account  the  veer  of  the  wind  with 
height,  especially  since  Paeschke  had  named  the  dimensionless  quantity  after  him. 
The  reply  by  Richardson  contains  a  number  of  interesting  comments  on  his  1920 

q 

paper  and  also  mention  of  his  further  explication  of  that  work  in  Richardson. 

The  Richardson  number  can  be  thought  of  as  the  ratio  of  the  destruction  of 
turbulence  by  stable  temperature  gradients  to  the  production  of  turbulence  by  wind 
shear.  Thus  turbulence  can  orcur  only  for  Ri  <  1.0.  This  criterion  is  now  widely 

Q 

used,  although  laboratory  tank  experiments  by  Thorpe  have  indicated  that  the  onset 
of  instability  does  not  occur  above  a  value  of  0.  25.  Also,  theoretical  studies  by 
Miles  and  Howard1^  suggest  that  Ric  -  0.  25  Is  a  necessary  but  not  sufficient  con¬ 
dition  for  the  onset  of  turbulence. 

Considering  the  flux  Richardson  number 


Ri 


Kh 

K~ 


Ri  . 


where  Kh  and  are  the  turbulent  coefficients  for  heat  and  momentum,  there  may 

be  some  critical  value  of  Hi.  <  1  where  instability  will  occur.  Perhaps,  as  pointed 

!  1  IC 

out  by  Vinnichenko  et  al,  under  certain  conditions  in  the  free  atmosphere,  either 
Kh  or  could  be  the  larger.  It  is  then  possible  that  instability  occurs  when 
RI  >  1.  In  this  report  it  is  assumed  that  Rlc  -  1  indicates  the  presence  of  turbulence. 
In  Figure  1,  the  frequency  of  occurrence  of  Ri  less  than  or  equal  to  0.  25,  0.40, 

0.  50,  0.  60,  0.  80,  and  1, 0  is  plotted  as  a  function  of  height  for  Brownsville,  Texas. 


6.  Paeschke.  W.  (1938)  Experimented  Untersuchongen  zum  Hauhigkelts  und 

Stabilltats  problem  in  der  bondennanen  Luftschicht,  Belt  range  zur  Physlk 
der  Atmosphare  (contributions  to  Atmospheric  Physics),  24:163-189. 

7.  Swinbank,  W.C.  ( 1952)  The  criterion  of  atmospheric  turbulence.  Quart.  J.  Roy. 

Meteor.  Soc.  Jjjj:420-425.  '  ’ 

8.  Richardson,  L.  F.  (1925)  Turbulence  and  vertical  temperature  difference  near 

trees,  Phil.  Mag.  S.  JC49,  81-89. 

9.  Thorpe,  S.  A.  (1968)  A  method  of  producing  a  shear  flow  in  a  stratified  fluid, 

J.  Fluid.  32:693-704. 

.  vvw 

10.  Miles.  J.W.,  and  Howard,  L.N.  (1964)  Note  on  a  heterogeneous  flow, 

J.  Fluid  Mech.  20:33  1-336. 

■  V4M 

11.  Vinnichenko,  N.  K.  ,  Pinos,  N.Z.,  Schimeter,  S.  M.  ,  and  Shur,  G.  N.  (1973) 

Turbulence  in  the  Free  Atmosphere,  Central  Aerologtcal  Observatory, 
Observatory,  Dolgoprudny,  USSR.  (Translated  from  Russian  by  John  A.  Dutton, 
Pennsylvania  State  University  for  the  Consultants  Bureau,  New  York.) 
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These  height  profiles  have  similar  structures;  however,  there  are  roughly  three 
times  more  occurrences  of  Ri  -  1.0  than  Ri  -  0.  25.  Data  from  several 
stations  were  examined  in  this  manner  and  yield  similar  results.  The  use  of 
Ric  S  1.0  provides  a  larger  data  base  for  statistical  analyses. 


Figure  3.  Frequency  of  Occurrences  of  Ri.  The  fitting  algorithm 
provides  a  smooth  fit  (solid  line)  through  the  original  east-west 
component  data  obtained  on  1C  February  1975  at  Brownsville,  Texas 


3.  DATA  BASE 

The  rawinsonde  data  used  in  this  study  have  been  provided  on  magnetic  tape  by 
the  National  Climatic  Center.  These  consist  of  twice  dally  "upper  air  observations" 
of  winds  and  temperatures  taken  at  0000  and  1200  hours  GMT.  Data  from  140 
stations  (mainly  from  the  United  States  and  South  America)  for  the  years  1971-1975 
(Figure  2,  Table  1)  have  been  reformatted  for  computer  processing  to  provide 
sample  statistics  for  the  Richardson  number  study. 

The  rawinsonde  system  was  implemented  to  provide  data  on  a  global  basis  at 
standard  millibar  levels  (mandatory  levels)  to  obtain  charts  and  maps  for  weather 
analysis.  A  calibrated  baroswitch  is  used  to  measure  pressure  levels  and  a 
hygristor  and  thermistor  are  used  for  the  measurement  of  relative  humidity  and 
temperature.  A  transmitted  signal  from  the  balloon  package  is  tracked  automatic¬ 
ally  to  determine  the  azimuth  and  elevation.  Altitude  and  winds  are  obtained  by 
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using  this  information  along  with  the  geopotential  height.  Hess  has  shown  that 
the  pressure  at  a  geopotential  height  can  be  determined  to  within  1  mb  by  using  a 
simple  barometric  element  of  the  type  used  in  the  rawinsonde  system. 

During  rawinsonde  measurements  a  person  in  the  data  recording  loop  records 
significant  temperature  levels  between  mandatory  levels  if  the  slope  of  the  traces 
on  the  recorder  record  changes  by  a  small  but  detectable  amount.  Significant 
temperature  levels  are  recorded  If  a  departure  from  the  trend  between  successive 
levels  represents  a  change  greater  than  ±  1°C  from  the  surface  up  to  100  mb 
(approximately  16  km)  and  greater  than  ±  2°C  up  to  the  level  of  termination.  The 
balloon  rise  rate  is  about  300  m  per  min  and  position  data  is  obtained  at  1-min 
intervals.  The  average  component  wind  speeds  are  obtained  from  the  position  data 
by  using  a  2-min  overlapping  interval  below  14  km  and  a  4-min  overlapping  interval 
at  higher  altitudes.  This  approach  results  in  wind  averaging  layers  of  approxi¬ 
mately  0.  6  km  below  14  km  and  1.  2  km  above  14  km.  A  compilation  of  standard 
error  estimates  for  meteorological  data  and  a  list  of  references  dealing  with  the 

subject  of  errors  in  meteorological  data  has  been  published  by  the  Meteorological 

13 

Group  Range  Commander's  Council.  In  the  present  report  it  is  assumed  that 
there  are  no  bias  errors  in  the  data.  Also  random  errors  are  minimized  by  group¬ 
ing  the  data  into  three-month  seasons  to  provide  the  largest  number  of  sots  of 
measurements.  This  number  ranges,  on  the  average,  from  180  at  the  2  km  level 
to  100  at  the  25  km  level.  The  number  of  measurements  decreases  as  altitude  in¬ 
creases  because  of  balloon  bursts. 

The  data  for  each  launch  were  read  and  verified  for  completeness  (Conversion 
to  standard  units  was  done  on  pressure,  altitude,  and  temperature).  To  obtain 
temperatures  and  gradients  of  temperatures  and  of  wind  components  required  to 
determine  the  Richardson  number  at  the  1-km  levels  used  in  this  study,  the  vertical 

profiles  of  wind  and  temperature  were  fit  using  a  standard  cubic  Hermlte  interpola- 

14  15 

tion  algorithm,  programmed  by  Tsipouras  and  Cormier,  This  algorithm  pro¬ 
vides  a  smooth  fit  passing  through  all  the  data  points  and  allows  the  computerized 
determination  of  the  fitting  function  and  its  derivative  at  any  point  while  avoiding 
problems  associated  with  discontinuities  encountered  in  linear  fits  (that  is,  a  first 

12.  Hess,  S.  L.  (1959)  Introduction  to  Theoretical  Meteorology,  Henry  Holt  and 

Company,  New  York,  N.V. 

13.  Meteorological  Group  Range  Commander's  Council,  1977:  Meteorological  Data 

Error  Estimates  (Document  110-77,  Secretarial  Range  Commander's 
Council,  White  Sands  Missile  Range,  New  Mexico  88002, 

14.  Forsythe,  G, ,  Malcolm,  M.  (1977)  computer  Methods  for  Mathematical 

Computations.  Prentice  Hall,  Inc.,  Englewood  Cliffs,  New  Jersey. 

15.  Tsipouras,  P. ,  and  Cormier,  R.  F.  (1973)  Hermlte  Interpolation  Algorithm 

for  Constructing  Reasonable  Analytic  Curves  Through  Discrete  Data  Points, 
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order  hold).  In  Figure  3,  the  original  east-west  wind  speed  component  data  from  a 
typical  set  of  measurements  are  overlaid  (continuous  curve)  with  the  output  of  the 
fitting  algorithm.  It  should  be  noted  that  there  is  substantial  resolution  in  the 
original  data,  usually  over  one  data  point  per  km,  implying  that  the  fits  should  be 
accurate.  Also,  the  1-km  intervals  chosen  are  consistent  with  the  original  data 
rate. 

In  Figure  4,  the  straight  lines  connect  at  points  which  represent  the  first  differ 
ences  in  the  original  wind  component  data  of  Figure  3,  at  altitudes  taken  at  the  mid¬ 
points  between  successive  data  points.  The  dots  in  Figure  4  are  the  output  of  the 
wind  shear  algorithm  of  the  fitting  function.  This  indicates  that  the  derivative  of 
the  fittmg  function  is  a  good  approximation  to  the  average  gradient  as  obtained  from 
the  first  differences  of  the  wind  component  data.  The  vertical  marks  at  the  bottom 
of  Figure  4  represent  altitude  levels  at  which  the  original  data  have  been  obtained. 


To  determine  the  average  resolution  of  the  data,  the  average  distance  between 
each  point  was  calculated  and  compiled  over  a  5-yr  seasonal  average.  Results  show; 

Altitude  Bin 

2-7.  5  km 
7.5-13.5  km 
13.  5-20.  5  km 
20.  5-25  km 

This  demonstrates  that  data  collection  was  consistent  (average  distance  be¬ 
tween  readings  falls  within  one  standard  deviation  of  their  mean)  and  Interpolation 
to  whole  kilometers  was  appropriate.  Variables  computed  for  1-kra  Intervals  from 
2  to  25  km  include  wind  shear,  density,  temperature  gradient.  Brunt -Vaisalii 
frequency,  turbulent  intensity,  rate  of  dissipation  of  kinetic  energy,  turbulent 
diffuslvity,  and  the  Richardson  number. 

Data  were  grouped  seasonally  to  assure  statistical  stability.  Seasons  were 
separated  using  the  solstice  and  equinox  dates.  There  were  from  100  to  1"0 
measurements  at  any  given  altitude  between  2  km  and  25  km. 

The  Richardson  number  is  computed  using  the  tempe rati. res  and  wind  compo¬ 
nent  values  along  with  their  derlvatlvt  i.  Percent  occurrences  of  Ri's  less  than  or 
equal  to  unity  were  computed  for  each  'eason  and  year  at  each  altituoj.  Values 
lesa  than  or  equal  to  1  are  used  as  ar  t  ••  ator  of  the  onset  tui  ence.  For  each 
station  there  are  980  data  readings  (7  years  X  4  seasons  per  year  -  35  altitude  bina). 

A  representative  set  of  21  stations  (Table  2)  were  chosen  from  the  -  stations 
In  the  continental  United  States  (Figure  5).  They  were  chosen  for  a  wide  range  of 
latitudes  and  longitudes  spanning  the  country.  Also,  sites  close  to  each  other  were 
included  to  measure  correlation  between  sites.  Topography  was  another  considera¬ 
tion.  It  was  desirable  to  include  stations  v  'th  varying  physical  configurations  to 
get  a  general  representation  of  the  country.  Years  1971-1975  were  chosen  since 
the  data  were  the  most  recent  at  the  time  and  most  of  the  sites  were  reporting  on  a 
regular  basis  through  this  period.  Also  stations  were  chosen  on  the  basis  that  they 
were  using  similar  and  up-to-date  equipment. 

The  number  of  samples  at  each  kilometer  bin  was  recorded  for  the  21  stations. 

A  minimum  of  100  samples  was  required  for  statistical  accuracy.  Table  3  lists 
the  highest  altitude  bln  that  attains  a  minimum  of  100  samples.  Over  85%  of  the 
groups  have  100  samples  at  25  km.  Due  to  the  geographic  formations  surrounding 
some  stations,  balloon  readings  began  at  2  km.  The  altitude  range  2  to  25  km  was 
used  for  analysis. 
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Table  2.  Stations  Used  In  Multiple  Regression  Analysis 


Station  Location 

Latitude 

Longitude 

Brownsville,  TX 

25°  54' 

97’  26' 

Chatham,  MA 

41*  40' 

69’  58’ 

Dayton,  OH 

39°  52' 

84’  07' 

Denver,  CO 

39°  46' 

104’  53' 

Flint,  MI 

42°  58' 

83“  44' 

Glasgow,  MT 

48°  13’ 

106’  37’ 

Green  Bay,  WI 

44°  29' 

CO 

o 

o 

CO 

CO 

Greenboro,  NC 

36°  05' 

79’  57' 

International  Falls,  MN 

48°  34' 

93“  23' 

Medford,  OR 

42*  22' 

122“  52' 

Miami,  FL 

25"  48' 

80’  16’ 

North  Platte,  NE 

41“  08' 

100“  41' 

Portland,  ME 

43“  39' 

70“  or 

Salem,  OR 

44“  55' 

123“  or 

Saulte  Ste.  Marie,  Ml 

46“  28' 

84“  22' 

Spokane,  WA 

47“  28’ 

117“  32' 

Topeka,  WA 

39“  04' 

95’  37’ 

Washington,  DC 

38“  59' 

77“  28' 

Waycr088,  GA 

31“  15' 

82“  24' 

Winslow.  AZ 

35’  01' 

110“  44' 

Great  Falls.  MT 

47’  29' 

111“  21' 

Between  seasons,  there  were  more  fluctuations  than  within  seasons.  For 
example,  height  profilt  s  of  turbulence  (Figures  7  through  27)  show  how  turbulence 
differs  seasonally.  Winter  consistently  has  the  highest  Richardson  number  range 
(that  is,  the  highest  turbulence  range)  and  the  greatest  fluctuations.  Conversely, 
3ummer  has  the  lowest  Richardson  number  range  (that  is,  lowest  turbulence  range) 
Thus,  the  data  available  has  been  reduced  from  980  data  points  to  a  more  con¬ 
sistent  480  points,  (24  altitude  bins  x  4  seasons  X  5  years). 
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Figure  5.  Location  of  the  81  Available  Rawinsonde  Stations  in  the  Continental 
United  States.  The  star  marks  those  stations  used  in  the  regression  analysis 

4.  regression  analysis 

Multiple  regression  techniques  were  used  to  investigate  the  patterns  of  tur- 

1 6 

bulence  over  the  continental  United  States.  Using  this  technique,  a  quantifiable 
relationship  was  established  between  the  location  variables  (latitude  and  longitude), 
altitude,  seasonal  and  yearly  variation,  and  the  dependent  variable  percent  Ri  less 
than  or  equal  to  one.  By  way  of  the  regression  coefficients  and  t  statist)"':  (ratio 
of  regression  coefficients  to  their  standard  errors),  we  can  quantify  the  relation¬ 
ship  of  turbulence  to  our  explanatory  variable  set. 

The  patterns  in  the  percent  occurrences  of  Ri  ^  1.  0  as  a  turbulence  indicator 
show  three  peaks  or  maximums;  one  below  7  km,  one  between  8  to  13  km,  and  one 
in  the  14  to  19  kin  range.  This  is  apparent  from  the  plot  of  Brownsville,  Texas 
(Figura  6).  in  order  to  account  for  these  peaks  in  a  regression  analysis  without 
having  to  use  higher  order  polynomials  in  altitude,  it  was  necessary  to  partition 
the  data  set  into  four  altitude  bins  and  analyse  each  of  the  regions  separately: 

2-7  km,  8-13  km,  14-19  km,  and  20-25  km. 

16.  Draper,  N. ,  and  Smith,  H.  (1966)  Applied  Regression  Analysis,  John  Wiley 
and  Sons,  Inc,,  New  York,  N.  Y. 
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Thus,  multiple  regression  was  the  analysis  used  to  quantify  turbulence  for  the 
21  stations  selected  in  the  continental  United  States.  The  altitude  region  2-25  km 
is  the  region  of  interest;  data  collection  is  consistent  over  this  range.  The  region 
is  divided  into  four  groups-.  2-7  km,  8-13  km,  14-19  km,  and  20-25  km  to  capture 
peaks  appearing  in  each  bin. 


Table  3.  Highest  Altitude  With  Greater  Than  100  Samples 


W 

SP 

SU 

F 

5'-year 

Average 

Brownsville 

25 

28 

27 

27 

26.75  km 

Chatham 

23 

26 

26 

25 

25 

Dayton 

24 

28 

28 

26 

26.5 

Denver 

26 

28 

28 

26 

27 

Flint 

23 

27 

25 

26 

25.  25 

Glasgow 

28 

28 

26 

24 

26,5 

Great  Falls 

24 

28 

28 

24 

26 

Green  Bay 

24 

21 

28 

26 

24.75 

Greensboro 

18 

26 

28 

26 

24.5 

International  Falls 

26 

28 

28 

26 

27 

Medford 

26 

28 

29 

27 

27.  5 

Miami 

27 

28 

26 

26 

26.  75 

North  Platte 

24 

28 

28 

26 

26.  5 

Portland 

18 

27 

28 

26 

24.75 

Salem 

26 

28 

28 

26 

27 

Sault  Ste.  Marie 

26 

28 

28 

28 

27.  5 

Spokane 

26 

28 

28 

26 

27 

Topeka 

24 

28 

27 

24 

25.75 

Washington.  D.C. 

26 

28 

28 

27 

27.  25 

Waycross 

26 

28 

28 

28 

27.  5 

Winslow 

25 

26 

26 

26 

25.75 

Average 

24.5 

27.  2 

27.4 

26.  0 

26.  3 

Figure  6.  The  Troposphere  and  Stratosphere  is  Divided  Into  Four 
Regions  for  Separate  Regression  Analysis.  The  Brownsville,  Texas 
occurrences  have  peaks  which  are  typical  of  the  peaks  in  these 
regions  generally  found  in  these  data 


In  symbols  a  multiple  regression  model  is  written  as  follows: 

*"  Vsjx2*-  0kxk.< 

where 

V  =  dependent  variable  (that  is.  To  Ri  5  1). 

=  ith  independent  variable  (four  seasons:  Winter,  Summer, 
Spring,  Fall;  (altitude)*,  i=  1,2,3:  (latitude)^ 

(longitudc)J  j  =  1,  2;  5-yearly  indicators  for  1971,  1972, 
1973,  1974,  1975), 

/IF  r  ith  regression  coefficient, 

t  *  random  error  term, 
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Yearly  variation  was  quantified  by  use  of  five  indicators  (zero-one)  variables 
(197  1,  1972,  1973,  1974  and  1975).  For  example: 


1971  = 


1  if  reading  taken  in  1971 


0  otherwise 


The  season  indicators  (winter,  spring,  summer,  and  fall)  are  defined  in  a 
similar  fashion. 


Winter  = 


1  if  reading  taken  in  winter 


0  otherwise 


Forthis  study,  a  stepwise  regression  analysis  is  performed.  Any  variable  is 
a  candidate  for  entry  into  the  model.  At  the  first  step,  the  variable  most  highly 
correlated  with  the  dependent  variable  is  entered  into  the  model.  Next,  the  partial 
F  statistics  associated  with  each  of  the  remaining  variables  are  computed.  These 
are  based  on  the  regression  equation  containing  that  variable  and  the  variable  ini¬ 
tially  selected.  The  variable  with  the  highest  partial  F  statistic  is  next  entered. 

This  process  continues  until  none  of  the  remaining  variables  have  a  significant 
partial  F  (that  is,  at  the  0.  05  level  of  Significance).  A  variable  that  has  entered 
may  be  removed  later  if  it  becomes  nonsignificant  with  other  variables  in  the  model 
(that  is,  if  the  partial  F  statistic  becomes  nonsignificant). 

Two  data  sets  were  developed: 

(1)  Four  reasons,  altitude1*  <k  =  1,2,3).  latitude-*  (j  =  1.2), 
longitude-*  (j  =  1,  2),  years  (1971-1975), 

(2)  Set  (1)  without  years. 

Since  data  is  ccllected  over  a  five-year  span,  set  (2)  w«»a  developed  to  examine 

o 

the  significance  of  years  in  a  regression  model.  If  the  R  (percent  variation  In 
distribution  of  percent  occurrences  of  Ri  -  1  explained  by  independent  variables) 
using  set  (1)  is  not  significantly  larger  than  the  corresponding  values  for  set  (2), 
we  can  Ignore  the  yearly  variation. 

The  amounts  of  variation  in  turbulence  are  increased  by  including  cross-product 
terms  of  altitude,  latitude  and  longitude.  The  increase  in  percent  variation  ac¬ 
counted  for  ranges  from  3-14%.  Examination  of  residuals  (actual  percent  occur¬ 
rences  of  Richardson  number  £  1  minuB  predicted  percent  from  the  regression 
including  cross  product)  indicated  that  improvements  in  the  model  could  be  made 
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in  the  8-13  km  range.  This  was  accomplished  by  the  inclusion  of  two  new  vari¬ 
ables  quantifying  season  and  year.  The  percent  variation  explained  by  the  inde¬ 
pendent  variables  increased  from  53%  to  60%  simply  by  adding  the  two  terms. 


5.  CHARACTERISTICS  OF  THE  VERTICAL  PROFILES  OF 
OCCURRENCES  OF  Ri  <  1 

Pronounced  seasonal  and  latitudinal  variations  are  evident  in  the  height  pro¬ 
files  of  percent  occurrences  of  Ri  s  1  shown  in  Figures  7  through  27.  The  general 
features  in  the  height  profiles  include  a  sharp  cutoff  in  the  occurrences  of  turbulence 
in  the  proximity  of  the  tropopause  level  and  the  varying  degree  of  turbulence  near  or 
just  above  the  planetary  boundary  layer  (above  2  km).  The  latter  feature  depends, 
of  course,  on  the  roughness  of  terrain  at  the  earth's  surface. 

At  latitudes  below  approximately  40°  N  latitude,  the  contrast  between  summer 
and  winter  occurrences  i9  evident.  At  t'  lowest  latitudes,  near  25°  N,  the  summer 
occurrences  are  as  low  as  5%  while  winter  approaches  40%.  The  spring  and  fall 
values  are  generally  about  equal  to  one  another  and  lie  between  those  of  winter  and 
summer.  As  latitude  Increases  northward,  the  number  of  occurrences  for  all 
seasons  converges  to  approximately  the  same  values.  Another  main  characteristic 
in  the  height  profiles  of  occurrences  in  turbulence  is  a  major  peak  in  occurrence 
near  the  10  km  altitude  level.  This  peak  occurs  in  all  seasons  but  is  much  more 
pronounced  in  the  winter  occurrences  at  low  latitudes  (40%)  and  is  lowest  in 
summer  at  low  latitudes  (10%). 

For  stations  near  mountain  ranges,  the  percent  occurrences  increase  above 
the  boundary  layer.  The  degree  of  turbulence  above  the  boundary  layer  appears 
to  be  related  to  the  height  of  the  mountain  ranges  in  the  proximity  of  the  data  re¬ 
cording  stations.  Height  profiles  obtained  from  data  at  these  stations  are  charac¬ 
teristically  quite  different  from  those  obtained  from  data  at  stations  located  in 
relatively  flat  regions.  Effects  from  the  mountains  appear  to  dominate  the  height 
profiles  of  occurrence  of  turbulence  to  the  degree  that  seasonal  features  as  well  as 
the  10  km  peak  are  no  longer  evident.  This  effect  i.i  apparent  in  the  profiles  of 
Glasgow,  Montana  (Figure  12),  Great  Falls,  Montana  (Figure  13)  and  Denver,  Colo¬ 
rado  (Figure  10)  where  the  highest  mountains  within  several  tens  of  miles  of  the 
station  location  are,  respectively,  3  K  ft,  7  K  ft  and  9  K  ft.  The  altitude  of  tur¬ 
bulence  occurrence  increases  with  height  of  surrounding  mountains. 
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6.  CONCLUSIONS 


Rawlnsonde  wind  and  temperature  data  have  been  used  to  determine  Richard¬ 
son's  numbers  for  1  km  levels  between  2  km  and  25  km.  A  critical  value  of  the 
gradient  Richardson  number  of  one  has  been  used  as  an  indicator  of  the  likelihood  of 
occurrence  of  turbulence.  The  patterns  in  the  occurrence  of  turbulence  have  been 
Investigated  for  the  seasons  and  years  from  197  1-1975.  A  regression  analysis  :■  - 
performed  on  data  from  21  rawlnsonde  stations  located  within  the  continental 
United  States  .indicated  that: 

(1)  Yearly  variations  in  the  occurrences  of  Ri  -  1  were  found  to 
be  small;  these  contributed,  at  most,  less  than  1%  of  the 
total  variation. 

(2)  The  amounts  of  variation  in  turbulence  explained  by  knowing 
latitude,  longitude,  altitude  and  season  included  38%  in  the 
2-7  km  bin,  60%  in  the  8-13  bin.  42%  in  the  14-19  km  bin 
and  11%  in  the  20-25  km  bin. 

(3)  The  patterns  of  occurrence  of  turbulence  based  on  a  critical 
Richardson  number  are  stable  and  reproducible  from  year  to  year. 

(4)  Examination  of  the  vortical  profiles  of  occurrences  of  Ri  ^  1 
reveal  that  there  are  pronounced  latitudinal  and  seasonal 
variations.  For  the  lower  latitude  data  {—  25 0  N),  the  summer 
occurrences  are  as  low  as  5%  while  winter  occurrences  approa 
40%.  Spring  and  fall  occurrences  are  about  equal  and  lie 
between  winter  and  summer  extremes.  Aa  latitude  increases 
northward,  the  number  of  occurrences  for  all  seasons 
converges  to  approximate. y  the  same  values. 

(5)  A  peak  m  occurrence  near  the  10  km  altitude  level  is  found 
In  the  height  profiles  from  low  to  mid-latitudes.  The  height 
profiles  of  occurrences  of  Ri  -  1  are  markedly  different  for 
stations  in  tne  vicinity  of  mountain  ranges. 
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PERCENT  OCCURENCE 


ore  7.  Seasonal  Height  Profiles  of  Occurrences  at  Ri  *=  1  for  Brownsville.  Texas  Station  Used  in 
Regression  Analysis 


rtftCCMT  OCCURENCE 


8.  Seasonal  Height  Profiles  of  Occurrences  at  Ri  £  1  for  Chatham,  Mass.  Station 
sion  Analysis 
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Figure  IS.  Seasonal  Height  Profiles  of  Occurrences  at  Ri  -  I  for  International  Falls,  MN  Station  Used 
in  the  Regression  Analysis 
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Figure  17.  Seasonal  Height  Profiles  of  Occurrences  at  Ri  ^  1  for  Miami.  Fla.  Station 
Regression  Analysis 
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